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Abstract 

We investigate how the background evolution affects the curvature perturbations 
generated by the curvaton, assuming a curvaton potential that may deviate slightly 
from the quadratic one, and parameterizing the background fluid density as p oc a~ a , 
where a is the scale factor, and a depends on the background fluid. It turns out 
that the more there is deviation from the quadratic case, the more pronounced is the 
dependence of the curvature perturbation on a. We also show that the background 
can have a significant effect on the nonlinearity parameters /nl and #nl- As an 
example, if at the onset of the curvaton oscillation there is a dimension 6 contribution 
to the potential at 5 % level and the energy fraction of the curvaton to the total one at 
the time of its decay is at 1%, we find variations A/nl ~ O(10) and A#nl ~ 0(1O 4 ) 
between matter and radiation dominated backgrounds. Moreover, we demonstrate 
that there is a relation between /nl and <7nl that can be used to probe the form of 
the curvaton potential and the equation of state of the background fluid. 



1 Introduction 



Although quantum fluctuations of the inflaton are often taken to be responsible for the 
origin of density perturbations, other mechanisms such as the curvaton [1-3], where fluc- 
tuations of a scalar field other than the inflaton generate primordial perturbations, have 
also attracted much attention recently. In particular, in the light of the recent result from 
WMAP5 which suggests that primordial non-Gaussianity may be large [4, 5]jj, the cur- 
vaton mechanism may be attractive since a large primordial non-Gaussianity can be gen- 



erated in this scenario [6-20 ^3, whereas simplest inflation models predict a non-linearity 
parameter /nl that is of the order of the slow-roll parameters (or 0(1) at most), and hence 
practically imply a Gaussian perturbation. 

Usually the curvaton potential is assumed to be quadratic. However, there is no other 
reason for this except simplicity, and in fact in any realistic particle physics model the 
curvaton can be expected to have some self-interactions. Thus, deviations from the exact 
quadratic form are also worth investigating. They have been discussed in [8,15-17,19-22], 
where it has been pointed out that non-quadratic contributions to the potential can modify 
the resultant curvaton perturbations in a significant manner. In particular, the prediction 
for the non-linearity parameters /nl and #nl can change considerably as compared to the 
quadratic case. 

In addition, there is yet another assumption which is tacitly adopted in the curvaton 
literature: the background evolution of the universe is determined by radiation. With this 
assumption, the curvaton starts to oscillate during a radiation-dominated (RD) epoch. 
If the curvaton decays before dominating the Universe, radiation is always the dominant 
component and controls the background evolution of the Universe. However, it is also 
possible that after inflation, the inflaton is oscillating around the minimum of the potential 
for a while and that the curvaton begins to oscillate during such epoch. In this case, the 
background evolution is different from the case of radiation and is determined by a matter- 



like component if the inflaton potential is approximatively quadratics 3 I. After inflation 
there could also exist a possibility of a kination-dominated phase where the kinetic term 
of some scalar field can dominate the energy density of the universe. Such fluid has a stiff 
equation of state with w = 1 while its energy density decreases as p oc a~ 6 . 

In this paper, we investigate how the background evolution of the universe affects the 
curvature perturbation generated from the curvaton. We first show that when the curvaton 
potential has a quadratic form, the background evolution has little effect on the curvature 
perturbation. However, when the curvaton potential includes a non-quadratic term, the 
nonlinear evolution of the curvaton field is affected much by the background, and the 



# x The degree of non-Gaussianity is usually characterized by the lowest order non-linearity parameter 
/nl- The current constraint on /nl is —9 < /nl < HI hi Ref. [4] or —4 < /nl < 80 at 95% C.L. in 
Ref. [5]. Although purely Gaussian fluctuations with /nl = are allowed, the central value is away from 
zero. 

# 2 Other models, for example, such as the modulated reheating scenario [23, 24] are also known to 
generate large non-Gaussian fluctuations [25-28]. 
# 3 In Ref. [18], this kind of situation is also included in the analysis for the quadratic curvaton potential. 
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resultant curvature perturbations can be significantly modified from the usual RD case. We 
investigate this issue by assuming a general background fluid and a potential that slightly 
deviates from a quadratic form, and derive the dependence of the curvature perturbation 
and/or the non-linearity parameters such as /nl and ^nl on different background fluids. 

The structure of the paper is as follows. In the next section, we summarize the formal- 
ism and give definitions of the quantities required for the subsequent analysis. In section [3], 
we discuss how the background evolution affects the quantities such as power spectrum 
and non-linearity parameters, using the appropriate formulas for arbitrary background 
fluids. The final section is devoted to a summary and a discussion of the results. 

2 Formalism and Definitions 

Let us begin by summarizing the formalism and the definitions of the various quantities 
required for the subsequent discussion. Here we consider a potential of the curvaton field 
a which, in addition to the usual quadratic term, also includes a non-renormalizable term: 



where m a is the mass of the curvaton, and A is a constant. For the purpose of this 
paper, it is enough to investigate the case of a slight deviation from the purely quadratic 
form. Thus in the following we assume that the quadratic term always dominates over the 
non-quadratic term (for a general discussion of the ramifications of non-quadratic terms 
in curvaton models, see [21,22]). We characterize the relative contribution of the non- 
quadratic term at the time when the curvaton is still in a slowly-rolling regime by the 
parameter s, defined as 



To investigate the curvature perturbation ( generated by the curvaton field, we adopt 
the 6N formalism [29-32] and calculate ( up to the third order as 



Once we obtain ( up to the third order, the power spectrum P^, bispectum B^, and 
trispectrum are given by 






(2) 




(3) 




(4) 



(C&CfeCfc,) = (^fBdh, k 2 , h)5(h + h + k 3 ). 



(5) 
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where and can be written as 



B c (h, h, k 3 ) = ~/nl (P ( (h)P c (k 2 ) + P c (h)P c (k 3 ) + P ( (h)P ( (h)) , (7) 

T f (fci, fc 2 , fc 3 , fc 4 ) = r NL (P f (fci 3 )P f (fc 3 )P f (fc 4 ) + 11 perms.) 
54 

+-S NL (P^P^P^h) + 3 perms.) . (8) 

Here /nl, tnl and ^nl are non-linearity parameters often used in the literature. Note that, 
in our case, r NL is related to / NL by 
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Thus, alternatively, we may assume the following expansion as the definition of /nl and 

C = Ci + |/nlCi + ^nlCi 3 + ■ ■ • , (10) 

where Ci denotes the curvature perturbation at linear order. 

Let us now derive the formulae for the power spectrum, /nl and #nl for a general 
background equation of state. As we mentioned in the Introduction, usually one assumes 
that the curvaton oscillates in a radiation background evolution. However, it is obvious 
that the inflaton can be expected to oscillate around the minimum of its (quadratic) 
potential for some time before it decays. For a weakly coupled inflaton field, the duration 
of this epoch could be fairly long. During that time the inflaton behaves like matter, and 
it is possible that the curvaton oscillations begin during this epoch. 

Eventually the inflaton will decay into radiation; this can well happen before the curva- 
ton oscillations dominate the energy budget of the universe, or before the curvaton decay, 
whichever happens first. Hence the universe becomes radiation-dominated so that the 
background evolution for oscillating curvaton is controlled first by a matter component, 
followed by radiation domination. Having this kind of situation in mind, let us assume 
there is a transition in the background from one fluid to another one at time t tT . We then 
parametrize the energy density of the background fluid before the transition time by 

Peg oc a~ a , (11) 

whereas after the transition, for t > t tr , we have 

Peg oca"' 3 . (12) 

For the case of oscillating inflaton background, the transition epoch t tT corresponds to the 
time when the inflaton decays. In this case, for t < t tT , a = 3 while for t > t tT , /3 — 4. 
Another example is a universe that is first kination-dominated so that the kinetic energy 
of a scalar field dominates the energy density of the universe. Then, after some time, 
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the universe becomes radiation-dominated. In this case, a = 6 for t < t tv while after the 
transition time (3 — 4. 

With this parametrization, we obtain the curvature perturbation at the linear order as 



Ci 



2cr ocs 
3(T sc 



R5a 



where 



R = r dcc (1 - kr tI ) + kr tr . 



(13) 



(14) 



The prime denotes the derivative with respect to er*. Here r dec and r tr roughly correspond 
to the fraction of energy density of the curvaton to the total energy density at the time 
of the curvaton decay and the transition of the background, respectively. Their precise 
definitions are 



3p 



apse + 3p CT 
Furthermore, k is defined by 



3p CT 



tr 



PPBG + 3p CT 



dec 



1 - 



a 

Tr 



(15) 



(16) 



By calculating the curvature perturbation up to the third order, we find that the non- 
linearity parameters /nl and g^L read as 



f: 



NL 



3 

2R 
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9 ^ Cr osc Cr osc g "osc^osc _^ 



cr 



A3 



0", 



9<j n sr d-R / 3 cr osc c" osc 
2 "KcF 



2# 3 rfa n 



9<c ^ 
4# 3 daL. 



(18) 

Here the derivatives of R with respect to cr sc appear. They are given in terms of the 
derivatives of r dec and r tr with respect to cr osc as 



dR 

deface 



1 - kr tI )- h ft(l - r de 



dcr n 



d(To 



and 



d 2 R 
da 2 



kr 



, d 2 r de 



tr 



dr^ec dr dec ^d 2 r tr 
2k- ; h k[l — r de 



d&osc dc osc dd osc d(J osc 
The derivatives dr dcc /da osc and drt T /da osc are explicitly given by 

dr dec 2 



da os , 
dr tr 



3cr n 



d(J osc 3cr osc 



-r d ec(l - ^dcc)[3 + ((3 - 3)r dcc ](l - fcrtr), 
-r tr (l - r tr )[3 + (a - 3)r tr ] 



(19) 
(20) 

(21) 
(22) 
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The second derivatives of r dec and r tr with respect to <t osc can be derived by differentiating 
the above equations and they are calculated as 



d 2 r dcc 



2 dr, 



dec 



3cr ns( , da n 



+ (1 - 2r dcc ){3 + (/3 - 3)r dcc }(l - kr tI ) 



+^dcc(l - r dcc )(p - 3)(1 - kr tr )] - -^-^^-r dcc (l - r dcc ){3 + (/3 - 3)r dcc }, 



3cT n sr d(J n 



(23) 



d 2 r tv 



dr 



tr 



3(7, 



3c d<j osc 



+ (1 - 2r tr ){3 + (a - 3)r tr } + (a - 3)r tr (l - r tr ) 



(24) 



We can now write down /nl and #nl as functions of the parameters a, (3, r dec and r tr in 
an explicit way, although in general the expressions are very complicated. In fact, as far 
as we consider only the cases with a, (3 > 3, r tr should always be smaller than r dec and in 
most cases we may assume that r tr <C r dcc . In this case, R ~ r doc and 0, /nl and <?nl may 
approximately be written as 



q ' dec 
<J C"osc 


5a*, 
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(25) 
(26) 



Cosc^osc 
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9(3 + (3 2 )+3r 2 dcc ((3-3) 2 



(27) 



As one can notice from these expressions, when r dec is small, which is the interesting 
case because then non-Gaussianity can be large, the effects of the background are almost 
encoded in the changes of ct osc and its derivatives. The quantities such as cr osc , cr osc and 
ct ' sc are evaluated at the beginning of the curvaton oscillation, and thus the background 
evolution after the transition is irrelevant if we assume a nearly quadratic potential for 
the curvaton. Thus in this case, after its change the background evolution does not 
affect much the curvature perturbation but modifies (j and its non-linearity parameters 
at most of order 0(1). This can be seen directly from the above expressions: (3 affects the 
coefficients by Oil). 

If there is no change in the background evolution up to or until after the time at which 
the curvaton decays, the expressions for the curvature perturbation and non-linearity 
parameters are quite similar to the case with r tr <C r dec mentioned above. In this case, we 
can simply set k = and R = r dcc in the above equations. Then d, / NL and #nl can be 
found by replacing (3 with a in Eqs. (I2"5l) . (I2"6"|) and fT2T|) . 
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3 Background evolution and non-Gaussianity 



Let us now apply the formalism of the previous Section and discuss how the background 
evolution affects the amplitude of the curvature perturbation and its non-linearity param- 
eters. As we pointed out, the transition in the background does not affect the perturbation 
much if the transition occurs much before curvaton decay. Hence here we focus on the case 
of no transition in the background, i.e., a single fluid controls the background evolution 
from the time well before the curvaton begins to oscillate to well after the curvaton has 
decayed. 

When there is no transition in the background evolution we can adopt the formulae 
Eqs. f j25|) . ff26l) and ff271) by replacing (3 by a in the equations, as mentioned above. To 
appreciate the impact of background evolution on the curvature perturbation and its 
nonlinearity, we plot £1, /nl and #nl in Figs. [THE] In Fig. the values of /nl and #nl 
are shown as a function of n (see Eq. ((2]) for the definition of the dimension of the non- 
quadratic contribution to the potential) for the cases of a = 3, 4 and 6 which correspond 
to matter, radiation and kination background, respectively, in our notation. As seen from 
the figures, the larger the non-quadratic power n is, the more pronounced the effect of the 
background becomes. 

To observe the magnitude of the effect more quantitatively, the values of Ci; /nl and 
#nl relative to those for the case with a = 4 and 3 are plotted in Figs.[2Hl]for several values 
of n as a function of a. In these figures, we fix the relative strength of the non-quadratic 
part s and r& cc to s = 0.05 and r^ ec = 0.01. In Figs. [5H7]we also plot the same information 
but this time for several values of s, fixing n = 8 and rd ec - 

Interestingly, when the potential of the curvaton is close to a purely quadratic one, the 
background evolution does not affect the results much. However, whenever the potential 
deviates from the quadratic form, the background tends to suppress the curvature pertur- 
bation d the more non-quadratic the potential is. It should also be noted that when a 
increases, or the background fluid becomes more stiff, £i becomes suppressed. For fixed s 
and n, the stiffening of the background fluid drives /nl and <?nl to increasingly negative 
territory. A point worth stressing is that these modifications are not small but e.g. by 
changing of the radiation dominated background to matter dominated one induces shifts 
in the non-linearity parameters that in principle could be easily observable. 

Finally, let us comment on the relation between /nl and #nl- As pointed out in [15], 
when the potential of the curvaton has a purely quadratic form and r^ ec is small, /nl and 
#nl are related as 

fl'NL ^ -y/NL- (28) 

However, when the potential deviates from a quadratic form, "the consistency relation" 
between / NL and <?nl becomes 

_ ^ _e2 /^osc^osc | o^oscOoscA ( -i ^oscCoscA /„„-> 
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where r^cc *C 1 is assumed. Because of the nonlinear evolution of cr, which is encoded 
in a osc and its derivatives appearing on the right hand side of Eq. (|29|) . the nonlinearity 
parameters are now related as — <?nl oc f^ L . This is in sharp contrast to the quadratic 
case where — #nl oc /nl- Furthermore, since the coefficient of the relation depends on the 
nonlinear evolution of a, it can be affected by the form of the potential and the background 
evolution. In this respect, the comparison of /nl and — S'nl, if they ever are observecF 4 



can be very useful for probing the form of the curvaton potential and the equation of state 
of the background during the time when the curvaton fluctuations generate the curvature 
perturbation. 

To demonstrate this explicitly, in Fig. M we plot the value of #nl as a function of /nl 
for several values of a with n = 8 and s = 0.02 . For reference, we also plot the case 
of the pure quadratic potential with RD background (a = 4). When / NL is large, which 
corresponds to the case of r^ec *C 1, the above relations hold. As a consequence, as one can 
see in Fig. El there is a definite difference in the #nl — /nl relation for different backgrounds 
even when the non-quadratic term in the curvaton potential remains the same. Thus the 
trispectrum may also provide important information not only about the curvaton potential 
but also about the equation of state of the background fluid during curvaton oscillations. 
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Figure 1: Plots of /nl and #nl as a function of n for several values of a. The values of s 
and r dec are taken as s = 0.05 and r dec = 0.01. 



^ 4 Recently, a constraint on (?nl has been obtained for the case of a negligible /nl; the limit is —3.5 x 
10 5 < 5 NL < 8.2 x 10 5 [33]. 
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Figure 2: Plots of Ci = Ci/Ci Ua ra > which is normalized to Ci for the pure quadratic case, 
relative to that for the cases with a = 4 (left) and a = 3 (right). The values of s and r^c 
are taken as s = 0.05 and r^ DC = 0.01. 
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Figure 3: Plots of /nl relative to that for the cases with a = 4 (left) and a = 3 (right). 
The values of s and r^ec are taken as s = 0.05 and r^c = 0.01. 
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Figure 4: Plots of #nl relative to that for the cases with a = 4 (left) and a = 3 (right). 
The values of s and r^ec are taken as s = 0.05 and r dec = 0.01. 
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Figure 5: Plots of £i = Ci/Ci Ua ™ > which is normalized to £i for the pure quadratic case, 
relative to that for the cases with a = 4 (left) and a = 3 (right). The values of n and rd ec 
are taken as n = 8 and rd ec = 0.01. 
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Figure 6: Plots of /nl relative to that for the cases with a — 4 (left) and a = 3 (right). 
The values of n and r^ec are taken as n = 8 and rd cc = 0.01. 





Figure 7: Plots of #nl relative to that for the cases with a = 4 (left) and a = 3 (right). 
The values of n and r^c are taken as n = 8 and r^cc = 0.01. 
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Figure 8: A relation between /nl and #nl for the case with a = 3, 4 and 6. The values of 
s and n are taken as s = 0.02 and n = 8. For reference, a pure quadratic case with a = 4 
is also shown. 

4 Conclusion 

In this paper, we have investigated the effect of the background evolution on the curvaton 
non-Gaussianity, assuming a curvaton potential which slightly deviates from the quadratic 
one. Such a study is motivated by the possibility that after inflation, the infiaton keeps 
oscillating about its global minimum for a long time so that the curvaton could actually 
decay while the universe is still effectively matter dominated. More exotic temporary 
possibilities, such as a kination driven universe, could also be envisaged. Therefore we 
have considered an ideal background fluid with some generic equation of state, leading to 
a background evolution oc a~ a , with a a free parameter. 

It turns out that the changing of the background to radiation, or equivalently, a change 
in the value of a, that takes place during curvaton oscillations, has by itself little effect on 
the perturbation. What matters is the nature of the background evolution before radiation 
domination finally kicks in, and as we show, it can lead to significant and potentially 
observable consequences. The non-linearity parameters /nl and #nl, as well as the linear 
curvature perturbation £i, depend on the nature of the background fluid. We find that 
the dependence on the background fluid becomes more pronounced as the deviation of 
the curvaton potential from the quadratic one increases. Typically, when replacing one 
background fluid with another, one induces effects on /nl that are easily of the order 
of O(10) but could also be much larger, depending on the relative strength of the non- 
quadratic part of the potential. 

An interesting issue is the relation between /nl and <?nl- We have showed that the 
relation depends both on the form of the curvaton potential and the background evolution, 
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which we find a rather surprising result. Hence measuring both /nl and #nl, or both the 
bispectrum and the trispectrum, would yield information not only on the curvaton self- 
interactions but also on the equation of state of the background fluid. Since that is linked 
to dynamics in the inflaton sector, here arises a possibility of probing inflaton physics at 
the very end of inflation. 
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